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EQUILIBRUM STATES FOR A CLASS 
OF DUAL-SPIN SPACECRAFT* 

by 

Thomas W. Flatley 
Goddard Space Flight Center 


CHAPTER 1 
INTRODUCTION 


Since the dramatic and unexpected attitude instability experienced by the first United States 
artificial satellite (Reference 1), now a textbook example (Reference 2) of the possible effects of en- 
ergy dissipation on “torque-free rigid body” motion, a design rule of thumb for gyroscopically stabi- 
lized spacecraft has been that such a body must spin about its maximum moment of inertia axis. 
Recently, however, a control concept that promises to relieve this increasingly burdensome constraint, 
while keeping much of the basic simplicity of the simple spinners, has emerged. It is known as dual- 
spin stabilization and in its simplest form features a rotor mounted on the nominal spin axis of the 
spacecraft and some type of passive damper. 

The basic idea had its beginning in the largely intuitive work of Vernon Landon in 1962. Unfor- 
tunately, his early claims were greeted with great skepticism, and his first attempt to publish his results 
was rejected. The dual-spin principle reached the open literature in 1963 (Reference 3) and 1964 
(Reference 4) and has been the subject of intensive research and frequent publication over the second 
half of the 1960’s. In fact, a symposium (Reference 5) held in El Segundo, California, in 1967 was 
devoted entirely to this subject. 

Understandably, published analyses have focused almost exclusively on the stability of pure rota- 
tions about the nominal spin axis. A variety of approaches to this problem are found in the literature, 
most of which involve linearization of equations of motion (derived either heuristically or exactly for 
specific systems) and analysis by informal averaging, “energy sink” approaches, and the more rigorous 
Routhian, Floquet, and Lyapunov methods. In some of these papers, one of the bodies is assumed to 
be completely despun. 


*This report was submitted as a thesis in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Aerospace 
Engineering, The Catholic University of America, May 1970. 
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However, the nature of dual-spin instabilities has been largely ignored in the literature. The sub- 
ject has been mentioned, notably by Fang (References 6 and 7) and Cloutier (References 8 and 9), but 
mostly in connnection with possible small system principal-axis shifts and equilibrium states involving 
damper masses displaced from their nominal positions. Cloutier (Reference 1 0) also discusses possible 
instabilities in the motion of a specific type of damper whose behavior is governed by Mathieu-type 
equations. 

A general analysis of dual-spin instabilities is not available, however. The research reported here 
addresses this problem by analyzing the detailed nonlinear equations of motion for a particular but 
fairly general dual-spin system by a perturbation analysis in which only the damper masses are assumed 
small and which involves variations in a general precessional motion of the system with an arbitrarily 
large cone angle. 

Among the various papers listed as references at the end of this report, the work of Mingori (Ref- 
erence 1 1) is particularly pertinent to the research reported here. The particular system chosen for 
detailed analysis is basically that studied in his paper, the notation used is similar to his, and the results 
obtained are frequently compared to those presented in the referenced publication. 

The basic objective of this work is to develop a better understanding of the physical principles 
involved in dual-spin stabilization by seeking equilibrium states for a dual-spin system (analogous to 
the “flat-spin” or “tumble” condition for single-spin bodies) which represent final states for unstable 
dual-spin configurations. 

Terminology can sometimes be a problem in applied mechanics, particularly in the area of rota- 
tional dynamics. Imprecisely or inconsistently defined terms and colloquial expressions can often be 
confusing. A brief review of the general motion of a torque-free symmetric rigid body will serve to 
define the various terms as they are used herein. 

In Figure 1.1, the 3-axis is the axis of symmetry for the cylindrical body shown and will be re- 
ferred to as the nominal spin axis. That component of the angular velocity vector parallel to this axis 
will be called the spin rate\ terms such as spin component and spin moment of inertia will refer to the 
nominal spin axis. The word transverse will imply a direction normal to the nominal spin axis (e.g., 
transverse angular velocity, transverse moment of inertia, and so forth). 

The angle 6 between the inertially fixed angular momentum vector H and the nominal spin axis is 
the nutation angle , or cone angle. In the strict sense, nutation involves variations in this angle; but, in 
accordance with popular usage and for lack of a better word, the rate at which the transverse angular 
velocity vector rotates in a body-fixed coordinate system will be referred to as the nutation rate. 

The motion of the nominal spin axis around the angular momentum vector is precession , and the 
corresponding rate is the precession rate. One source of confusion is the fact that this phenomenon is 
observable only when there is a nutation angle, i.e., if the angular momentum vector and the nominal 
spin axis are not coincident. In this case, the motion is popularly referred to by various authors as 
precession , nutation , coning , or wobble . 
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H For a symmetric body, the total angular ve- 

locity vector and the nominal spin axis are co- 
planar with the angular momentum vector and 
rotate about it at the precession rate. 

With a rigid body, all of the above-mentioned 
quantities are constants of the motion. It is now 
well known, however, that inevitable energy dissi- 
pation in real systems (i.e., quasi-rigid bodies) can 
drastically change the characteristics of the mo- 
tion. The effect is critically dependent on the 
moment-of-inertia ratio q = C/A, where C and A 
are the spin and transverse moments of inertia, 
respectively. For q > 1, any coning is damped 
out, and an eventual pure rotation about the nom- 
inal spin axis is attained. With q < 1, the nutation 
angle increases with time toward a stable equilib- 
rium value of 90°, reaching what is commonly 

called a flat spin or tumble. In any case, the final state is a minimum-energy state for the system con- 
sistent with the amount of angular momentum present. 
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CHAPTER 2 


SYSTEM DEFINITION AND EQUATIONS OF MOTION 


The particular dual-spin system to be analyzed consists of two primary bodies with a common 
axis of rotation and a motor that maintains an arbitrary but constant relative spin rate a. For conven- 
ience, these bodies will be referred to herein as a spacecraft and a wheel , although the relative size, 
weight, and moments of inertia of the two parts are unrestricted. Each is basically symmetric with 
respect to its common axis, which will be called the nominal spin axis. Schematically, they will be 
represented as bodies of revolution (Figure 2.1). The spacecraft has a mass M , a spin moment of iner- 
tia / 3 , and a transverse moment of inertia I y Corresponding properties of the wheel are M\ I r v and I' v 
The centers of mass of the two bodies are a fixed distance / apart. 

Attached to each of the primary bodies is a spring-mass-dashpot damper with a single degree of 
freedom parallel to the nominal spin axis. The spacecraft damper is at a distance a from the spin axis 
in a plane containing the center of mass. It has a total mass m b , an active mass m, a spring constant k , 
a damping coefficient c, and a displacement coordinate z. The wheel damper is similarly specified by 
the constants a\ m f h , m\ k\ and c' and a displacement coordinate z\ 

To provide static and dynamic balance of the system 
when the active masses of the dampers are at rest (i.e., 
z = z' = 0), six balance weights, three each of masses m b and 
m by are added as shown in Figure 2.2. This figure also shows 
a set of transverse principal axes for each of the bodies, arbi- 
trarily chosen such that the /- and /'-axes pass through the 
dampers. 

Figure 2.3 shows a complete description of the system 
to be analyzed. 

The equations of motion for this system are presented, 
without derivation, by Mingori in Reference 1 1. They are 
rederived in Appendix I and agree in every detail with those 
in the literature. They form a seventh-order, nonlinear, 
nonautonomous set, inconveniently coupled even in the 
highest derivatives. 

In essence, the first three equations (2.1 through 2.3) 
govern the basic attitude behavior of the system and are 



Figure 2.1— Dual-spin system. 
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Figure 2.2— Balance weights and transverse principal axes. 


analogous to the classical Euler equations of motion for a rigid body. Equations 2.4 and 2.5 are 
second-order equations which primarily determine the motion of the active masses of the dampers. 

In spite of the complex nature of the system of equations, we are assured a priori of the existence 
of a first integral of the motion corresponding to the total angular momentum. In principle, this fact 
could be used to reduce the order of the system by one, but of more interest here is the fact that the 
lack of external torques guarantees that the total angular momentum vector will be a constant of the 
motion. We will take advantage of both the constant magnitude of this vector and the inertially fixed 
reference direction it provides. 

[A +M r $ 2 + mz 2 - 2m(X + l 2 )z + m'z' 2 - 2 m'(? ~ l^z'jaj x - ( maz + m'a'z' cos \jj) co 3 + m'a' sin \ p z' 
= (A - C)gj> 2 co 3 - J' 3 o(j0 2 - 2M t ^co { + M T ^ 2 co 2 co 3 + 2 m{£ + / 2 )(ico 1 - zgj> 2 cj 3 ) 

+ mz( 2fcu 1 - 2zto 1 + zco 2 co 3 + aco l oo 2 ) + 2 m f {$ ~ l^iz'ooy - z'co 2 co 3 ) 

+ ra'z'(2f C 0 j - 2z t co 1 + z'co 2 co 3 + a' cosi// tOjO^) - m'a'z' sin V/ [(co 3 + cr) 2 - co 2 ] , (2.1) 
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3 , 3 ' 



[A + mz 2 - 2 m(f + / 2 )z + m f z' 2 - 2m'(f “ /^z'Jd^ ~ m ' a ' z ' sin \i/ cb 3 - maz - mV cos \p z' 

= (C - /l)a; 1 cc 3 + / 3 aco 1 - 27kf r ffa; 2 - Af r f 2 0 J 1 co 3 + 2m(X + l 2 )(zco 2 + zco 1 co 3 ) 

+ raz(2fco 2 - 2zco 2 - zco 1 co 3 ) + maz{ to 2 - co 2 ) + 2ra'(f - / 1 )(z'cc 1 co 3 + z'co 2 ) 

+ m'z'( 2fa> 2 - 2z'gj 2 - z'a^co^ + ra'z'[( co 3 + a) 2 - co 2 ]a' cos V/ - raV sini// z / co 1 cj 2 , (2.2) 
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- {maz + m'a'z' cos \p)co l ~ m'a' sin \p z'tb 2 + Cco 3 


= ma(2zco 1 - zco 2 co 3 ) + m'a f sini//(2z'c o 2 + z'co 1 co 3 ) + mV cos \p(2z r co l - z'co 2 co 3 ) , (2.3) 

-maco 2 +m(l - p)z - m'pz' --maco i co 3 + m{co\ + co 2 )[z(l - p) - 1 2 - p'z'] - cz-kz , (2.4) 

and 

m'a' sin \p gj x - m'a' cosi p co 2 - mp'z + m'(l - p')z' 

= - m'a'[sin i// co 2 (co 3 + 2a) + cos \ p cu^cOj + 2a)] 

+ m'(co 2 + co 2 )[z'( 1 _ p') + l x - pz] - c'z' - k'z f . (2.5) 

These equations are written in terms of the following useful combinations of constants and vari- 
ables in order to reduce the number of terms: 

M t = M + AT + 4 m b + 4m f b , 


m 



f = pz + P'z' 

/(Af + 4m^) 

= ’ 

/(M' + 4mp 

/ 2 = ‘ ’ 

a 

A = /j + /j + 2m b a 2 + 2 m b a' 2 + (Af + 4m^)(Af' + , 


C = / 3 + / 3 + 4 m b a 2 + 4m f b a f 2 , 

and 

J’ 3 = 1^ + 4m' b a’ 2 . 

The constant Af r is the total mass of the system, ,4 is the combined transverse moment of inertia 
of the system about the center of mass with the damper masses at rest, and C is the combined spin 
moment of inertia, which includes / 3 (the spin moment of inertia of the wheel, its damper, and its 
balance weights). 
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CHAPTER 3 


PERTURBATION ANALYSIS, PART I: 

UNPERTURBED MOTION OF THE SYSTEM 
AND FIRST-ORDER DAMPER MOTION 


An approximate solution of the equations of motion will be carried out with the assumption that 
the effects of the dampers are a small perturbation on the general “rigid body” motion of the system. 
The “zeroth-order” equations (which represent the unperturbed state, or the reference motion, of the 
system) are easily obtained if the masses m and m' are allowed to vanish in Equations 2.1, 2.2, and 2.3. 
The resulting set of equations is as follows: 


ACo x + (C - x4)gj 2 oo 3 + J\ooj 2 = 0 , 

(3.1) 

A6o 2 - (C - A)co 1 oj 3 - J 2 o(x> l = 0 , 

(3.2) 

Cgj 3 = 0 . 

(3.3) 


Equation 3.3 has the obvious solution co 3 = constant, and we will label that constant £2. Substi- 
tuting this solution into Equations 3. 1 and 3.2 reduces them to a pair of constant-coefficient linear 
equations with an elementary solution: 


and 

where 


co 3 + Xco 2 = 0 
co 2 - Xcjj = 0 , 


AX = (C-A)i2 +J' 3 o . 

Thus, co 3 and co 2 are sinusoidal functions with a constant amplitude co and frequency X. 
loss of generality, we can select a reference time when co 2 = 0 and write these solutions as 

co 3 = co cos Xt 
and 


co 2 = co sin Xt , 


Without 


where co is the amplitude of co 3 at t = 0. 
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A. Characteristics of Reference Motion 

The equations above show that the unperturbed motion of the dual-spin system being studied is 
very similar to the general precessional motion of the torque-free, symmetric, rigid body in classical 
mechanics. Figure 3. 1 is a representation of the behavior of the angular velocity vector in the spacecraft- 
fixed coordinate system. 

The 3-component (or the spacecraft spin rate) is constant with a magnitude of SI, and the trans- 
verse component is of fixed magnitude co but rotates around the 3-axis (or nominal spin axis) at a 
rate X. 

Figure 3.2 shows the motion of the angular momentum vector in the same reference frame, the 
spin and transverse components of angular momentum, and the cone angle d. 

The quantities co, SI, X, H, and 6 are constants of the unperturbed motion of the system. The 
actual motion of the nominal spin axis (Figure 3.3) is a steady precession on a cone of half-angle 6 
around a fixed angular momentum vector. This would be the motion seen by an outside observer, and 
the precession rate is given by \p 1 = H/A, obviously another constant of the motion. 


Qualitatively, this motion is exactly that of the classical, torque-free, symmetric, rigid body, but 
there are quantitative differences due to the presence of the wheel. The nutation rate X and precession 
rate \jj x are given by 



Figure 3.1 -Angular velocity vector in Figure 3.2-Angular momentum vector in Figure 3.3— Actual motion of 

body coordinates. body coordinates. nominal spin axis. 
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When o=0, the system is reduced to exactly the classical rigid body with the familiar expressions 




d 


and 


B. First-Order Damper Motion 


*1 = 


CD. C 

2B — 

A cos 0 A 


D for 6 < l 


To find a first approximation for the motion of the active masses in the two dampers, we turn to 
Equations 2.4 and 2.5, retain all terms linear in the damper masses m and m', and substitute the 
zeroth-order angular velocities determined above. The resulting equations are 


and 


mi + c'z + (k- mu> 2 )z = maco(X - D) cos Xt - mca 2 l 2 
m'z' + c'z' + ( k ' - m'<jj 2 )z' = m'a'to(X - D - 2a) cos (a - X)t + m' cu 2 /j . 


(3.4) 

(3.5) 


Equations 3.4 and 3.5 are similar in form and are simply representations of classical forced, 
damped, harmonic oscillations. Their steady-state solutions can be written as 


z = Cj + C 2 sin Xt + C 3 cos Xt 


and 


z’ = C 4 + C 5 sin (a - X)t + C 6 cos (a - X)t , 
where the constants C. can easily be shown to be 

co 2 /., 


C i = 


GJ 


- 


2 _i ’ 

m 


a — \co(X - £1) 
m 


c 3 = 


(— - co 2 - X 2 V + (— x\ 

\m / \m / 

aco(X - "Co 2 - X 2 ) 


(--co 2 -X 2 ) +(-X) 
\m / \m / 

c o 2 L 


2 ’ 


c 4 k' 


K 2 

C 


m 
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and 


C 5 = 


a ,± 7 co(a-X)(X-«-2ff) 

m 


[5-^-(X-a)f + [f,(a-X)] 2 ' 


a'co(X - £2 - - co 2 - (a - X) 2 J 

— , - co 2 - (X - a) 2 ] 2 + K (a - X)1 2 
Lm J L m J 


Thus far, we have determined the unperturbed motion of the dual-spin system and, in effect, have 
assumed that this basic motion drives the two dampers. Having found the steady-state response of the 
damper masses to the precessional excitation, we now seek the effect of their oscillatory motions on 
the attitude motion of the system. 
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CHAPTER 4 


PERTURBATION ANALYSIS, PART II: 

ATTITUDE DISTURBANCE AND 
STABILITY CONSIDERATIONS 


Since the physical system under consideration is free of external torques, the total angular mo- 
mentum vector will be a constant of the motion. The magnitude of this vector, H = |H|, will be a first 
integral of the exact equations of motion The fixed direction of the angular momentum vector pro- 
vides a useful inertial reference line, and in consideration of the attitude behavior of the system, it is 
convenient to concentrate attention on the motion of the nominal spin axis relative to this fixed line. 

The cone angle 0 < 6 < 7r between the angular momentum vector and the nominal spin axis is a 
key parameter relative to the stability of possible motions. Knowledge of its behavior as a function of 
time is the usual objective in stability analyses for spin-stabilized spacecraft. 

Now, if the constant angular momentum vector is decomposed into components perpendicular 
and parallel to the nominal spin axis (the 3-axis in Figure 4.1), the magnitude of the latter vector will 
obviously be 


H Z =H cos0. (4.1) 

For the system under consideration, however, this component of the angular momentum can also 
be shown to be 


H 3 = Cco 3 + / 3 a - mazWy - cos \p + co 2 sin \p) . (4.2) 

Thus, 

H cos# = Cco 3 +/ 3 a- mcfzcoj - m'a'z'ico^ cosi p + co 2 sini p) . 

Differentiating both sides of this equation with respect to time and eliminating co 3 by use of 
Equation 2.3 yields an exact relationship: 

-H sin d d = ma(zcOy - z« 2 co 3 ) + m'a'z'( co 3 + aXojj sin \p - co 2 cos i p) 

-t-m'u'z'fcoj cos^ + co 2 sini//) . (4.3) 
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Figure 4.1— Decomposition of 
angular momentum vector. 


Equation 4.3 represents an exact equation of motion for the system 
cone angle 8. The right side of this equation contains the small parameters 
m and m'. If m = m' = 0, the cone angle is a constant of the motion (as 
seen also in Chapter 3). 

To determine the first-order perturbation of the cone angle, we sub- 
stitute the approximate expressions for the body rates and the damper 
motions found in Chapter 3, repeated here for convenience: 


and 


Wj = co cos Af , 
co 2 = co sin A t , 
co 3 = £2 , 

z = C x + C 2 sin A t + C 3 cos A t , 
z' = C 4 + C 5 sin (a - A )t + C 6 cos (a - A )t . 


Substituting first for the body rates, assuming, without loss of generality, that \p(0) = 0 we obtain 


-H sin 8 8 = maoo{z cos A t - zSl sin A 0 + m'a'co[z'(£2 + a) sin ( a - A )t + z' cos (a - A)t] . (4.4) 

Now, substituting for z and z' would produce sin, cos, sin 2 , and cos 2 terms with arguments A t and 
(a - A )t. Through use of trigonometric identities, the squared terms may be eliminated, and the result- 
ing equation would contain only constants and sinusoidal oscillatory terms. In the spirit of the method 
of averaging (Reference 12), then, the systematic change of 6 is considered to be controlled by the 
constant terms, with the periodic terms contributing only small oscillations of 8 about some average 
value. 


With only the constant terms retained, the secular behavior of the cone angle is found to be gov- 
erned by the equation 


C C 

-H sin 8 8 = macu(A - £2) + m'a'u)(2o - A + J2)~ 


(4.5) 


Substituting the calculated values of C 2 and C 5 from Chapter 3 and approximating the transverse 
momentum H sin 6 by the transverse angular momentum of the unperturbed motion Alo yields 


8 = - 


H_ sin 8 j u 2 cA(A f2) 2 

2y4 2 


a' 2 c'(o - A)(2a - A + tt) 2 


- co 2 - A 2 )' + aV r^-co 2 -(a-A)2| 2 + K(a-A)' : 

\\m / \m / m J L m' J 


(4.6) 
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Now, each term within the brackets above is either a fixed constant or a constant of the unper- 
turbed motion of the system. These latter constants can be expressed explicitly in terms of the cone 
angle 0 as follows: 


co = 


H sin 0 
A ’ 


and 


n = 


\ = 


H COS 0 - / 3 CT 


(C - A)H cos 0 + AJ'j a 
AC 


Equation 4.6 is thus of the form 0 = /(0), and for a first-order differential equation such as this, 
questions of stability are almost trivial. Equilibrium will obviously be possible whenever /(0) = 0, and 
the stability of any equilibrium point depends on the value of a derivative of /(0), usually the first, at 
that point. Figure 4.2 shows this in graphical form. 

The labels U, M, and S identify points of unstable, metastable, and stable equilibrium, respec- 
tively; the arrows show the direction of motion along the phase plane trajectory. Stable equilibrium at 
0 = 0 Q is seen to require /(0 O ) = 0 and /'(0 O ) < 0. 

Equation 4.6 is rather remarkable in that it represents a single first-order differential equation 
that determines the overall attitude behavior of a seventh-order, nonlinear, nonautonomous system 
with arbitrary initial conditions. Plotting the function /(0) for any particular dual-spin configuration 
would immediately reveal the stability characteristics of the system, indicate “zones of attraction” for 
particular stable equilibrium points, and show the “time constant” associated with cone-angle buildup 
or decay near equilibrium states. 

The time constant r is given by the inverse of the slope of /(0) at the equilibrium point, as shown 
below: 


With 


0 = m 



Figure 4.2— Equilibrium and stability if f(d) = 0. 
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and 


f(9 Q ) = 0, 


let 

0 = 0 Q + x . 

Then, 

0 = x 

and 

m = f(d Q +x)~f(e 0 )+xf'(e 0 ), 

so 

x~x/'(0 o ). 

Thus, 

f =/’«>„>*• 

and the solution can be written as 

x = x(Q)e^' (e ^ . 

Since such exponential behavior is conventionally written as 

x = x(0)e ^ T , 

we have 

1 

T '/'(0 O )' 


It is convenient to let 0 = F sin 0, where 

/ 


F = 


H 

2 A 2 


a 2 cX(X - £l) 2 


T2 + 


a'V(a-X)(2a-X + ft) 2 


l(^-co 2 -X 2 ) + (^X) g-cuMa-X) 2 ] +[~(a-X)]J 


In this problem, 0 < 0 < it. Equilibrium always exists at 0° and 1 80° and possibly exists at inter- 
mediate angles if F vanishes. It follows that the dual-spin system is stable— 


(1) at 0 = 0° if F(0) < 0. 

(2) at 0 = 180° if F(it) > 0. 


(3) at 0 = 0 O if F(0 O ) = 0 and F'(0 O ) < 0. 

An immediate example will clarify these concepts. Consider the special case where there is a 
damper on the spacecraft only. Then, 
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H 

a 2 c\(\ - ft) 2 

2 A 2 

km-" 2 - X2 ) 2+ (m X )j 


and since only the sign of F is significant, the positive constants, the term (X - ft) 2 in the numerator, 
and the entire denominator can be factored out. Thus, F“-X, so that 9 « -X sin 9. 

Now, 


(C - A)H cos 9 + AJ'^a 


Hence, X(0) is a biased cosine function. In the special case of a = 0 (i.e., no relative motion of the 
wheel) and C> A, the function 9(9) would be of the form shown by the solid line in Figure 4.3. 

Note the stability at 0° and 180°; since the wheel has been despun and the system reduced to an 
ordinary spinning spacecraft, this result is in agreement with the well-known fact that rotation about 
the nominal spin axis is stable if the moment-of-inertia ratio is “favorable” (i.e., if C> A). If the case 
C <A had been considered, the function 9(9) would be the mirror image of that shown by the solid 
line and there would be unstable equilibrium at 0° and 1 80° and stable equilibrium at 9= 90°. This 
latter point would represent a flat-spin condition, the famous minimum-energy state for spinning 
satellites with “unfavorable” moment-of-inertia ratios. 

The dashed lines in Figure 4.2 show the effect of the addition of a positive relative wheel rate to 
the system. The unstable equilibrium point is moved to the right monotonically as the spin is increased, 
eventually reaching 1 80° as a critical rate a Cf is attained. For a > a cr , only one stable equilibrium 
angle (0°) exists, and the 1 80° point becomes unstable. An initial condition of 9 - 1 80° would then 
result in an eventual complete inversion of the system; i.e., the spacecraft would literally turn upside 


0 
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down in reaching a stable equilibrium state, a can be calculated by setting X = 0 and cos 6 = - 1 in 
the equation for X above; i.e., 


( C-A)H 

~aJ[ 
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CHAPTER 5 


ANALYSIS OF RESULTS 


A. The Fundamental Result Restated 

The primary result obtained in the preceding chapters can be summarized briefly as follows: The 
condition 6 = d Q represents a state of stable equilibrium (or a stable steady-state cone angle) if /( 0 O ) = 0 
and /'( 0 O ) < 0 , where 


f{d) = sin 6 F(6) , 


F(d) = 


H 

2 A 2 


a 2 c\(\ - Cl) 2 


a’ 2 c'(o-\)(2o-\ + n ) 2 


C~ -co 2 ~ X 2 ^ + (^x) 7 - co 2 - (a - X) 2 1 + K(a-X)] 

\m ) \m / m J m J ; 


(5.1) 


£2 = - 


H cos d- Fa 


X = 


(C - A)H cos 6 + AJ\a 


AC 


and 


co = 


H sin 6 


The function F represented by Equation 5.1 is a function of 14 variables, 1 1 of which would be 
fixed constants for a given spacecraft configuration (the moments of inertia C, A, and / 3 , and the 
damper parameters m, k, c, a, m', k', c', and a'). The other independent variables are the angular mo- 
mentum H, the wheel rate a (constants of the motion in this analysis), and the cone angle 6. The vari- 
ables Cl, X, and co are not independent, and they are given explicitly above. 

The most important quality of F, however, is its sign, and the various terms involved are predom- 
inantly positive. The only terms possibly negative are X and ( 0 - X) in the numerators of the two 
fractions. In terms of the independent variables, these quantities are 


X = 


(C- A)H cos d + AJ' 3 o 
AC 


(5.2) 
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and 


(A-C)H cosd +A(C-J')a 

°' X = AC' ‘ (53) 

Also, in the special cases of damping on one body only, the sign of F is determined by ( 1 ) the 
sign of (-X) for spacecraft damper only, and (2) the sign of (a - X) for wheel damper only. 

In what follows, initial states of the system will be considered to be characterized by a pure rota- 
tion about the nominal spin axis, with an initial spacecraft spin rate of and a wheel rate of o. An 
S2 0 a-plane then contains points representing all possible initial states for a particular system. The 
angular momentum corresponding to a particular point in this plane will be given by H = |CT2 0 + J^o\, 
and since Cf 2 Q + J'^a = H cos 0 Q , an initial state for which CSl Q + J'^a > 0 will represent an initial cone 
angle of 0°. Likewise, if C£l Q + J' 3 a < 0, the initial cone angle will be 1 80°. The line C£2 Q + J'^a = 0 
thus divides the £2 Q a-plane into two parts, as shown in Figure 5.1. 

B. Special Case I : Damper on Spacecraft Only 

The stability of an arbitrary initial state, with damping on the spacecraft only, is determined 
according to the following conditions: 

( 1 ) 0° is stable if X 0 is positive. 

(2) 1 80° is stable if X 0 is negative. 


a 
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In both cases, X Q is calculated from Equation 5.2, with the appropriate values of H and cos d, as follows: 

(C - ,4)(O2 0 + J’ 3 o) + AJ\a 
X ° = AC 


C^A 

A 





The line X Q = 0, or a/n Q = (A- C)/J' v then divides the initial state plane into regions of X Q > 0 
and X Q < 0, as shown in Figure 5.2. 

Combining Figures 5.1 and 5.2 and recalling the stability criterion above yields Figure 5.3, a sta- 
bility diagram for the system with damping on the spacecraft only. Regions marked S and U are stable 
and unstable, respectively. 

A useful alternate form of this stability diagram can be constructed in terms of dimensionless 
parameters. Let q = C/A and q' = J f 3 /A. The stability boundaries above then become the lines 
q + q\o/^l 0 ) = 0 and q + <jr'(a/J2 0 ) = 1. Figure 5.4 shows the new stability diagram. The only values 
of q that are physically meaningful lie between the limits q ' and 2. 

Figure 5.4 verifies a result obtained by Mingori (Reference 1 1) through a Routhian analysis of 
linearized equations of motion representing this special case and pure spins about the nominal spin 
axis. 


Turning now to the nature of possible instabilities, we recall that stable equilibrium states can 
exist at intermediate cone angles (0 < < 7r) if F(6 0 ) = 0. Thus, in the special case at hand we are 

interested in the possibility of X vanishing, where, in terms of dimensionless parameters, 


X = 


, t .Hcos8 , 
(<7“ 1 ) 2 — + qo 


(5.4) 



Figure 5.2— The sign of Xq. 
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where 

"j= \qSl Q +q'o\ . 

Now, X is seen to be a biased cosine function of the cone angle 6, and by inspection of Equation 
5.4 the following may be deduced: 

(1) For a = 0, X will vanish at 8 = tt/2 and will have a slope at that point determined by the sign 
of (q~ 1). 
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(2) Introducing a 0 monotonically raises or lowers the curve, moving the point where X van- 
ishes to the right or left and eventually eliminating the possibility of X vanishing for a real 0. Critical 
values of o are determined from Equation 5.4 with X = 0 and cos 0 = ± 1. 

These critical values, o 1 and a 2 , are calculated as follows: 

(q - l)(q£l Q + q f G x ) + q , o 1 = 0 , 
or 


or 


(1 - q)(q£l 0 + q'o 2 ) + q f o 2 = 0 , 

g 2 _ g(l -g) 
q\q - 2 ) * 


(5.6) 


Since /(0) = F(d) sin 0, the function /(0) will take on one of the four characteristic shapes shown 
in Table 5.1, depending on the magnitude of the wheel rate o and the sign of (q - 1). In this table, sta- 
ble equilibrium states are marked by an arrow (/) and comments are added concerning conditions that 
produce the various possibilities. 

Class I behavior includes as a special case the “favorable” single-spin spacecraft, which is stable 
right side up (0 = 0°) or upside down (0 = 1 80°). Note that only in Class II can a stable intermediate 
cone angle exist. The “unfavorable” single-spin spacecraft is a special case of this type of behavior, 
where 0 Q = 90°. Equation 5.4 can be used to calculate the value of the stable intermediate angle in 
general, with the result being 


q’o 

cos 0 — ■ 

0 (\- q )\q£l 0 +q'o\- 


(5.7) 


For Class II, all pure spins about the nominal spin axis are unstable. Classes III and IV include 
the possibility of a heretofore unreported type of instability, a complete inversion of the system. Un- 
stable pure spins about the nominal spin axis will change to stable pure spins about the nominal spin 
axis, with the spacecraft turning upside down in the process. 

Table 5.1 shows that only stable initial states and inversions are possible if q > 1 or, equivalently, 
if C > A. A complete stability diagram for this case is shown in Figure 5.5. 

For q < 1, two types of initial-state instabilities are possible— inversion, and a steady-state inter- 
mediate cone angle. One obvious observation is that if o = 0, the latter possibility will occur and the 
steady cone angle will be 90°. In general, for q < 1, Table 5.1 shows that stable intermediate angles 
exist if 

<7(1 -q) a ^ 1 ~<7 
q'(q - 2 ) £ 2 0 q' 
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le behavior for f(Q). 
Stable 

Cone Angles 


0 °, 180 ° 


0 ° < 0 O < 180 ° 



Conditions 


q>\ 


q < 1 


q < 1, a < a 7 


C/ * Q ' . ° 

S > ' / 

0 


spacecraft only and C > A . 



The upper bound here is actually the line q + #'O/^ 0 ) = 1 , which has been found to separate sta- 
ble and unstable initial states. The lower bound represents a new line, 

a = gO -q) 

q\q-2Y 

which divides the previously shown unstable regions into two parts corresponding to the two possible 
types of instability. Intermediate cone-angle values can be calculated from Equation 5.7, and the loci 
of points representing fixed values of d 0 will be straight lines in the £2 0 <j, initial-state plane. Figure 5.6 
shows a complete stability diagram for the case of q < 1 . 

Presenting the above information in the alternate form [in oji 2 Q , <?-space] has disadvantages and 
compensating advantages. Many of the straight lines in Figure 5.6 become curved, and the distinction 
between supplementary cone angles (i.e., 6 and 7T - d) is lost, but the effect of the parameter q is more 



60 - 


90° _ STABLE 

INTERMEDIATE 


ANGLES ^ ^ 


120 — 


150 


INTERMEDIATE 

ANGLES 


180 



* 120 ° 


150° 


Figure 5.6— Complete stability diagram for damper on the spacecraft only and q < 1. 
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clearly seen. For definiteness of cone angle only, it is assumed in the construction of Figure 5.7 that 
£2 0 > 0. For clarity, only that portion of the stability diagram involving the two types of instability is 
shown. Refer to Figure 5.4 for the overall view. 

C. Special Case 1 1 : Damper on Wheel Only 

The stability of an arbitrary initial state, with damping on the wheel only, is determined accord- 
ing to the following conditions: 

( 1 ) 0° is stable if X Q - a is positive. 

(2) 1 80° is stable if \ Q - a is negative. 

In these conditions, X Q - a is calculated from Equation 5.2 with appropriate values of H and cos 9 as 
follows: 


Q 
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X 0 -a = 


(C - A)(Cn Q + J' 3 a ) + A(J' 3 - C)a 
AC ~~ 


C-v4 

A 


S2 0 + 


J' z -A 

A 


a . 


The line A 0 - a = 0, or 


a C - A 

f2 0 A-jy 


then divides the initial-state plane into regions of - a > 0 and A 0 - o < 0, as shown in Figure 5.8. 

Combining Figures 5.1 and 5.8 and recalling the stability criterion above yields a stability diagram 
for damping on the wheel only, shown in Figure 5.9. Regions marked S and U are stable and unstable, 
respectively. 

Again, a useful alternate form of this stability diagram can be constructed in terms of dimension- 
less parameters. With q = C/A and q ' - J\/A , the stability boundaries above become the straight lines 


and 




« + «'#- =1+ ^r 


This alternate stability diagram, Figure 5. 10, verifies a result obtained by Mingori (Reference 1 1) 
with a Routhian analysis of linearized equations of motion representing this special case and pure spins 
about the nominal spin axis. 


o 
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Turning now to the nature of possible instabilities, we recall that stable equilibrium states can 
exist at intermediate cone angles (0 < d Q < n) if F(d 0 ) = 0. Thus, in the special case at hand we are in- 
terested in the possibility of X - a vanishing, where, in terms of the dimensionless parameters, 
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(5.8) 


and 


(q~ 1)^J cos0 + (q 1 - q)o 

q 

J= |?£2 0 + q'a\ . 


An analysis of Equation 5.8 would closely parallel the earlier analysis of Equation 5.4. Each is a 
biased cosine function of the cone angle 0 which vanishes at 0 = 90° if a = 0; its slope is determined by 
the sign of q - 1, and the bias level is a monotonic function of a. However, critical values of a are now 
determined for X - o = 0 in Equation 5.8, with cos0 = ±1. 


or 


or 


(q - 1 )(q£l 0 + q'o x ) + ( q ' - q)o t = 0 , 

q l _ q - 1 
n 0 1 -q'- 

( 1 - q)(q£l 0 + q'a 2 ) + ( q ' - q)o 2 = 0 , 

°_2 qi}_ -<?) 

~ q + qq’~ 2 q' ' 


(5.9) 


(5.10) 


The discussion of the behavior of the function /(0) with damping on the spacecraft only is appro- 
priate for this special case as well. The four classes of behavior represented in Table 5. 1 are again pos- 
sible. The conditions producing these classes are different, however. Table 5.2 shows the appropriate 
inequalities involved. 

All comments in the previous section concerning the various classes of behavior apply here also, 
except that in Class II, the value of the stable intermediate cone angle is determined from Equation 5.8 
to be 


cos 0„ = 


(<?' -q)o 


0 (1 -q)\qSl Q + q'o\' 


(5.11) 


Table 5.2 shows that only stable initial states and inversions are possible if q > 1 . A complete 
stability diagram for this case is shown in Figure 5. 1 1. 

For q < 1, we again have two possible types of instability— inversion, and a steady-state interme- 
diate cone angle. The latter possibility (Class II of Table 5.2) occurs when 

1 q gO ~0) 

i - q' <i + qq'- 2 q' ' 

The lower bound here is seen to be the line which has previously been found to separate 
regions of stability and instability in the initial-state plane above. The upper bound is a new line which 
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Table 5.2— Four classes of possible behavior for f(6). 



*The terms a x and u 2 are given by Equations 5,9 and 5.10, respectively. 


Stable 

Cone Angles 


0°, 1 80° 


0° < e 0 < 1 80° 


0° only 


1 80° only 


Conditions 


Q > 1 

0 2 < O < (Jj* 


q<\ 

o l <o <o 2 


q> 1 , o <o 2 
or 

q < 1 , a < a 2 


q < 1, a > Oj 
or 

q > 1, a > a 2 


divides the unstable region in Figure 5.9 into two parts, corresponding to the two possible types of 
instability. 

Appropriate values of intermediate cone angles can be determined from Equation 5.1 1, and the 
loci of points representing constant values of will be straight lines in the J2 0 a-plane. Figure 5.12 
shows the complete stability diagram for the case of q < 1 . 

The above information can also be presented in an alternate form (in the a/f2 0 , g-space) as shown 
in Figure 5.13. Again, for definiteness of cone angle only, the assumption is made that 0. Note 

the singularity in the boundary separating inversions from the intermediate angles. It occurs when the 
denominator of the expression 
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Figure 5.1 1— Complete stability diagram with damper on the wheel only and 

q> 1. 


a q( 1 ~ q) 
q + qq' ~ 2q' 

vanishes, or when 

2 q' 


Figure 5.13 shows only that portion of the (o/S2 Q )q plane of most interest here. (Refer to Figure 
5.10 for the overall view.) 

D. General Case With Tuned Dampers 

For the special cases considered so far, specific damper parameters have not been important. We 
have been concerned only with the moments of inertia C, A, and J f v the angular velocities £2 and go, 
and the cone angle Q. The general case, however, involves eight additional variables which appear in 
the function F (Equation 5.1 )—m, k , c, and a and m\ k\ c\ and a f . However, the stability characteris- 
tics of a particular configuration can easily be determined by a numerical evaluation of the function 
F(d ) and an examination of the results in light of the basic stability criteria we have been using. 

In order to simplify the analysis of the general case, we make the following assumptions, which 
reduce the number of independent parameters: 
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K 

i 


q 



\ i 

>1 — — a/J2 0 


Figure 5.13— Portion of stability diagram in alternate form, with damping on wheel only and £2 0 > 0. 

and 


We have seen several times that the stability of initial states is determined according to the fol- 
lowing conditions: 

(1) 0° is stable if F(0)<O. 

(2) 1 80° is stable if F(tt) > 0. 

These requirements can be reduced to one, viz., an arbitrary initial state is stable if F(0 Q ) cos 0 Q < 0, 
where cos0 Q = sgn(<?£2 0 + q'o). Thus, for stability of an initial state, it is required that 



33 












or, equivalently, 


(q & o + q'o) 


■ (X^V 2 (2a-X 0 + n 0 )^ 

' ■' ■ (a-x 0 )|a-x/ 


X 0 IX 0 I 


o' 


> 0 , 



(5.13) 


Examination of Equation 5.13 reveals that there exist four possible situations, which are shown 
in Table 5.3. 

Mingori (Reference 1 1) studied linearized equations representing the special case being considered 
here, with five particular combinations of the parameters q , q\ and a/f2 Q . These equations had peri- 
odic coefficients, and he explored a 5<5'-space at various isolated points by a numerical, computerized 
Floquet analysis, producing a total of 1 25 data points, to which the results obtained here can now be 
compared. 


Table 5.3— Initial state stability, damping on both bodies. 




Requirement for Stability 


+ 

f> 

(h\( 2 JL-_hL +l ) 

\n 0 /v n o ) 

l 2 

+ 

(h. \(jl 





L\n 0 /\n 0 nj 

- 


— 


Always stable 



+ 


Never stable 




f< 

\n 0 A n„ 7 

2 



(h. _ A 





_\^ 0 J\n 0 sij_ 
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Table 5.4 shows the particular combinations of parameters examined by Mingori and the corre- 
sponding stability criteria deduced from the analytical results obtained here. 

Figure 5.14 shows comparative results for Cases 1, 2, and 3 of Table 5.4, in which an analytical 
stability boundary is predicted. Nearly total agreement occurs, with the exception of single points in 
Cases 1 and 2. (These particular points will be discussed later.) 

Results for Case 4 are not shown; however, they agree with Mingori’s results, which yielded insta- 
bility at each point tested. 

In Case 5, Mingori found instability at over half of the points tested (Figure 5. 1 5), a result appar- 
ently in gross disagreement with the stability criteria reported here. To assess the significance of this 
disagreement, we digress to consider the work of Mingori in some detail. 

Recall that this result comes from a Floquet analysis, which was necessary because the linearized 
equations of motion had periodic coefficients with damping on each of the bodies. However, Mingori 
also presented the results of a Routhian analysis for the special cases of damping on one body only, 
which can be used to investigate the stability of the particular configurations examined on the 5- and 
5'-axes, i.e., when 5 or 5' vanishes. Appendix II contains a detailed discussion of the behavior of the 
system on these axes and shows damper size limitations which arise when the dampers are tuned. 

Briefly, the Routhian analyses produce five inequalities that are necessary and sufficient for sta- 
bility in the special cases to which they are applicable. Two of these are trivially satisfied for physi- 
cally realistic systems; a third represents that criterion which has been shown to be in agreement with 
the results obtained here, i.e., the stability diagrams presented earlier for those special cases. 

The other two inequalities pertain to the sizes of the damper masses, and Table 5.5 shows the 
damper size limitations imposed by these inequalities for the special combinations of variables exam- 
ined by Mingori. The headings II- 1 5, 11-16, and so forth refer to equations in Appendix II. 

Of particular interest in this table are the very low limits on the size of the spacecraft damper in 
Case 5. Spacecraft damping has been found to be stabilizing in this instance, yet instability is indicated 
for 5 >0.008. 


Table 5.4— Mingori’s particular combinations of q , q\ and a/f2 0 . 


Case 

q 

q' 

o 

£ 2 0 

q+q 'ik 

ii 0 

V 

J°I Q 

I' 

Requirements 
for Stability 

i 

0.8 

0.2 

4 

1.6 

0.6 

3.4 

5/5' > 13.7 

2 

0.8 

0.2 

-6 

-0.4 

-1.4 

-4.6 

5/5' > 1.48 

3 

0.8 

0.2 

-2 

0.4 

-0.6 

-1.4 

5/5' <0.414 

4 

0.8 

0.2 

0.5 

0.9 

-0.1 

0.6 

Never stable 

5 

1.5 

0.2 

-1.5 

1.2 

0.2 

- 1.7 

Always stable 


*\ 0 /n 0 =<7+?Wtto)-i. 
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Case 1 



6 

Case 2 



5 

Case 3 


KEY 

O - MINGORI'S STABLE POINTS 
■ - MINGORI'S UNSTABLE POINTS 


ANALYTICAL STABILITY BOUNDARY 


Figure 5.14— Comparison of results. Cases 1, 2, and 3. 


The inequalities involved in these low limits (derived in Appendix II) are repeated here for con- 
venience: 

\r 5 S W) + V >0 

and „ , , 

In an effort to find the physical significance of these inequalities, the appropriate linearized equa- 
tions in Reference 1 1 were examined. Terms extraneous to our purpose here were removed for clarity, 
and the following set of equations was found to possess possible pertinent types of instabilities: 
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0.1 Cr 

0.05 O' 


co x + A Q to 2 = 0 , 


(5.14) 



<5 

Figure 5.1 5— Mingori's Case 5 results. 

Aq > S£2q’ or equivalently, Aq - 5£2g < 0. This 
Routhian inequality above. 


"2-V>l =! j n S Z . <5.15) 

and 

mz + kz = ma(cj 2 - figCOj) . (5.16) 

Substituting Equation 5.15 into Equation 5.16 
produces 

s + (£- sa o)z = “(\- a o)“i- < 5 - 17 > 

Note that the “effective spring constant” in 
Equation 5.17 involves 5, and in order for this to be 
positive, it is necessary that 5 < 

If the damper is tuned to its excitation fre- 
quency, i.e., if k/m = Aq, this inequality becomes 
latter inequality is seen to be similar to the first 


If a fourth-order equation in z is formed, equivalent to the set of Equations 5.14 through 5.16, it 
is found that the linear term has the coefficient \*(k/m) - 5A Q £2j*. Requiring this “spring constant” to 
be positive implies kjm > 5fi^/A 0 . If the damper is again tuned to its excitation frequency, this re- 
quirement becomes Aj* > 5£2^/A Q , or equivalently, Aq - SAq^j* > 0. This inequality resembles the 
other Routhian inequality listed above. 


Physically, then, these inequalities involve effective spring constants which are characteristic of 
the coupled equations. In the perturbation analysis reported here, the equations of motion are un- 
coupled, and these types of instabilities could not be detected. 


Since the inequalities involved depend on the size of the dampers, it seems reasonable to consider 
them as representative of second-order effects beyond the scope of this perturbation analysis. The 


Table 5.5— Special case limits on damper sizes. 


Case 

Q 

q' 

o 

J2 Q 

5 

m 

11.15 

ax 

11.16 

m 

11.17 

ax 

11.18 

i 

0.8 

0.2 

4 

0.225 

0.178 

1.45 

No limit 

2 

0.8 

0.2 

-6 

No limit 

1.57 

10.6 

No limit 

3 

0.8 

0.2 

-2 

0.900 

No limit 

No limit 

1.56 

4 

0.8 

0.2 

0.5 

0.011 

No limit 

0.267 

No limit 

5 

1.5 

0.2 

-1.5 

0.033 

0.008 

4.83 

1.03 
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above inequalities show that tuning the dampers to low excitation frequencies can greatly enhance 
these effects and severely limit the size of the damper. This is the problem encountered here. 

If Table 5.5 is used as an indicator of the range of applicability of our first-order perturbation 
analysis, the following conclusions may be drawn: 

(1) None of the wheel damper size limitations ( 8 ' max ) are of any consequence. They are all re- 
mote from the range of primary interest (0.001 to 0.1), except for Case 4 and Equation 11.17, which is 
not important since there is first-order instability in that case. 

(2) The limitations on 6 in Cases 2 and 3 are probably remote enough from 0. 1 to be inconse- 
quential, and the low limit (0.01 1) in Case 4 is not significant because of the first-order instability in 
that case. 

(3) Both of the 6-limitations in Case 1 lie just beyond the range investigated by Mingori, so that 
the second-order effects may have destabilized the point of disagreement mentioned earlier. 

(4) In Case 5, the line 6 = 0.008 cuts through the middle of the stability diagram shown in Figure 
5. 14. The stability at all points to the left of this line is in agreement with the stability criterion devel- 
oped here. At one point, the stabilizing wheel damping was efficient enough to stabilize the system 
even though 6 exceeded 0.008. 

The point of disagreement in Case 2 remains unexplained. A situation similar to that of Case 1 
exists here, but the second-order instability boundary is further away (at 5 = 1.57). An alternate ex- 
planation of the disagreement at this point may lie in a sign error in a very small term in one of the 
matrices worked with by Mingori.* This error appears in two published versions of his work (Refer- 
ences 5 and 1 1), but may or may not be significant.** 

E. Transformation of the Stability Diagram 

The stability criteria developed here for the case of tuned dampers with 0. 1 critical damping on 
each body can be used to trace the transformation of the stability diagram as energy dissipation capa- 
bility is gradually shifted from one of the bodies to the other. The extreme cases would be the dia- 
grams shown in Figures 5.4 and 5.10 for damping on the spacecraft only and wheel only, respectively. 

Recall that the stability criteria for initial states, given by Equation 5.13, involved the signs of the 
quantities 


q q iT’ 


*The term W x (4, 3) should be negative. 

**In a recent personal communication, Mingori acknowledged the presence of this error in his published equations and in his computer 
program. He also brought to the writer’s attention, however, a paper by Lindh and Likins [“Infinite Determinant Methods for Stabil- 
ity Analysis of Periodic-Coefficient Differential Equations,” AIAA J 8(4), p. 680, 1970] in which results in total agreement with his 
are presented. He added that the same error is present in the equations in this paper but not in the computer program which produced 
the results. 
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I 



and 

a 


Figures 5.16, 5.17, and 5.18 show which regions of the ( o/£l 0 )q -, initial-state plane represent positive 
and negative areas for these quantities. 

Inspection of Table 5.3 reveals that the relative amount of damping on the two bodies is unim- 
portant if X 0 /^ 0 and o/£l Q - X 0 /f2 Q have opposite signs, i.e., all initial states are either stable or unsta- 
ble. Figure 5.19 shows the regions of the initial-state plane to which this situation applies. 

In other areas of this plane, the sign of the function F(0 Q ) given by Equation 5.12 will be impor- 
tant, and stability boundaries (or places where F changes sign) will be characterized by the vanishing 
of the bracketed quantity in Equation 5.13. Since we are now interested strictly in cases where A 0 /f2 0 
and a/£2 Q - \/^ 0 have the same sign, we can thus concentrate on satisfying the relationship 


-2. - . ,JL - i 

“o n o 


0 G 

, + , 


q 
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Figure 5.19— Special regions of the initial-state plane. 


8 _ V 2a x o + 2 

8' ~ [(X 0 - fi 0 )(a - X 0 )_ 

With j3 = (5/.5') 1/2 > solutions are sought for the equation 


A 0 ( 2 a - A 0 + n Q ) = ±p(X 0 - n Q )(a - A 0 ) . 

With the substitution A 0 = (q - 1)£2 0 + q'a, this equation becomes a quadratic function in a/Cl 0 
that can be solved analytically for lines along which F vanishes. Appendix III shows a derivation of 
the general quadratic and the particular form it takes when q' is arbitrarily selected to be 0.2. This 
particular case is used for illustrative purposes. The function F, referred to several times below, is re- 
peated here for convenience: 


F°c -5 


- sy 2 

x 0 ix 0 i 


+ 5' 


(2a ~ A 0 + S2 0 ) 2 
(a-A 0 )|a-A 0 | ' 
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Study of the transformation of the stability diagram is begun with an examination of the behav- 
ior of the function F with damping on the wheel only. The initial stability diagram is shown in Figure 
5.10 and is repeated in Figure 5.20. Two important lines corresponding to X Q = 0 and 2a - X Q + £2 Q = 0 
are also shown. Near the first of these lines, a small amount of spacecraft damping will strongly affect 
the function F; near the second, the stabilizing effect of the wheel damper is weak (i.e., F » 0). 

Now, the addition of a small amount of damping on the spacecraft splits both the stable and un- 
stable regions along the dashed lines in Figure 5.20. The result is shown qualitatively in Figure 5.21. 

The spacecraft damping has changed the sign of F in two regions, stabilizing one and destabilizing 
the other. A further increase in spacecraft damping causes these new regions to grow, completely 
eliminating the leftmost stable area above when j3 = 1. Adding still more damping causes the leftmost 
F= 0 line to sweep to the left, leaving stable area behind, as shown in Figure 5.22. This moving line 
eventually stabilizes the remainder of the left half-plane and then begins to approach from infinity on 
the right, leaving stable area behind again. As the spacecraft damping becomes more and more domi- 
nant, the stability diagram takes on the appearance of Figure 5.23. 

The final state, as the transformation is completed, is shown in Figure 5.24. Lines corre- 
sponding to A 0 - f2 Q = 0 and a - X Q = 0 have been added to this figure, indicating places where the 

stabilizing effect of the spacecraft damper is 
weak (i.e., F » 0) and where the effect of a small 
amount of damping on the wheel is strong, re- 
spectively. 

2 a - X 0 + £2 0 = 0 

Figure 5.25 shows the progressive steps of 
the transformation for the special case of q f = 0.2. 
In order to show the complete transformation on 
one fairly simple figure, many lines have been 
omitted and the figure appears somewhat abstract, 
° /il o but the information shown is quantitatively ac- 
Figure 5.20— Initial stability diagram. curate. 
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Figure 5.25— Transformation of stability diagram 




CHAPTER 6 


POSSIBLE GENERALIZATIONS 


In Chapter 3, it was said that, in effect, the unperturbed motion of the system was determined by 
letting the damper masses vanish; that motion was allowed to drive the dampers, and then their steady- 
state response to the excitation was found. In Chapter 4 it was shown that, for small damper masses, 
the spin component of the system angular momentum can be written as 

H cos 6 » Coo 3 + J^o , 

so that 

-H sin 6 6 « Cch 3 . 

Thus, cone-angle behavior is intimately related to spacecraft spin-rate changes. 

The damper motion mentioned above generates reaction torques on the spacecraft; in particular, 
it generates a spin torque that can change the spacecraft spin rate. From the sign of the average spin 
component of the reaction torque, the sign of th 3 , and consequently the sign of 0, can be determined. 

In order to assess the significance of the specific damper orientation chosen in this analysis and in 
the hope of expanding the generality of the results obtained, we consider here the effects of a general- 
ized spring-mass-dashpot damper (arbitrarily positioned with respect to the center of mass and arbi- 
trarily oriented relative to the nominal spin axis) on the motion of a symmetric dual-spin system, using 
the above-mentioned heuristic approach to the dynamics. 

The basic system consists of a symmetric spacecraft and a symmetric wheel driven at a constant 
relative angular velocity. Figure 6. 1 shows this system and gives several pertinent parameter definitions. 

A. Unperturbed Motion of the System 

The characteristics of the unperturbed motion of the system are given by the solution of the well- 
known Euler equations: 


-j-H + co T x H = 0 , 
dt T 

where 

H = Ao) + (Cfl + 
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UNIT ♦ 
VECTOR k 



is the total angular momentum vector and all vectors are measured in a spacecraft-fixed reference 
frame with origin at the center of mass and one axis along the nominal spin axis. 


With the appropriate substitutions and differentiation, 


or 


and 


Ago + cm + (w + m) X [A 03 + (02 + /' a)k] = 0 , 
Ago + (02 + J'^a - ASI)go x k = 0 


02 = 0 . 


Thus 


where 


£2 = constant and co = Xk x co , 


X = 


( C-A)Sl + J' 3 o 

A 


( 6 . 1 ) 


which is a constant. 
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If co = (cop co 2 , 0), Equation 6.1 has the elementary solution co = (co cosXt, co sinXt, 0), where co 
is a constant of the motion and we have chosen co 2 (0) = 0 without loss of generality. 

B. Location and Orientation of the Damper 

Let r Q and e specify the location and orientation, respectively, of a small, single-axis, spring-mass- 
dashpot damper with a spring constant k and a damping coefficient c. If p locates the active mass of 
the damper relative to its undeflected position, the actual position of the damper mass is given by 
r = r Q + p, where p = p^e. Without loss of generality, we can set one of the components of r 0 equal to 
zero. We will then be working with the vectors 


and 




P 




C. Damper Equation of Motion 

The motion of the active mass of the damper will be given by the solution of the vector equation 
ma = F, where m is the mass, a is the total acceleration, and F is the total force. In terms of the above 
vectors, this becomes 

F = m[r + 2co r x r + co T x r + co r x (oj t x r)] 

= m[p + 2 oj t xp + cbxr + cjj.x (oj^ x r)] 

= m[p + 2c o T x p + X(k x cu) x r + u> t (oj t • r) - r(w|.)] 
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= m[p + 2 o3 t x p + Xou(r • k) - Xk(r • co) + <x> T (co T • r) - r(co|,)] . 


( 6 . 2 ) 


The equation for the axial motion of the mass is obtained from the scalar product of each side of 
Equation 6.2 and the vector e, with the right side becoming ~kp x - cp v Thus, 

m[p 1 + X(e • «)(r • k) - Xe 3 (r • co) + (e • a> r )( cj t • r) - co|.(e • r)] + cp x + kp x = 0 . (6.3) 

D. Spin Component of Reaction T orque on Spacecraft 

A reaction force - F acts at r, so the resulting torque on the spacecraft is T = F x r. The system is 
assumed to be balanced such that this torque vanishes when the damper is caged, i.e., when p 3 = pj = 0. 
If F q is the corresponding force, it is required that T„ = F fl x r Q = 0. Of primary concern is the spin 
component of the torque T - T Q , or 


r 5 = k-(Fxr-F 0 xr 0 ) 


= k-[Fx Pl -r 0 x(F-F 0 )], 


where F is given by Equation 6.2 and 

F - F 0 = m[p + 2oj t x p + Xco(p * k) - Xk(p • co) + io T (io T • p) - poo^. ] . 

Appendix IV gives the derivation of the final result: 

T s = m{- e 2 x Q p l + Pj [2c 3 (r • 0 J r ) - 2£2(p 3 + r Q • e)] 

- p 1 (e l 0J 2 - e 2 c OjIfXIZg + p 3 e 3 ) + 03 T • r] - p 1 x 0 co 2 (Xc 3 + c o T ■ e)| . (6.4) 

E. Special Case I: Damper Parallel to Nominal Spin Axis 

We now let e l = e 2 = 0 and e 3 = 1 . The resulting equations are 

Pi + ^Pi + (^‘" c ° 2 ) p i = n)x 0 oj cos\t + z Q co 2 

and (6.5) 

T s = mx 0 co[2p 1 cos Xt - p 3 (X + S2) sinXf] . 

The first of these equations is identical to Equation 3.4, produced by the perturbation analysis. 
The variables P 1 ,x Q , and y 0 are the same quantities labelled z, a, and - 1 2 previously. This equation 
can be solved in closed form, and if the steady-state solution is substituted into the second of Equa- 
tions 6.5, an average value of T s can be determined in a way analogous to that which produced an 
average value of 6 in Chapter 4. 


I 
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The equations above can be manipulated to represent the motion and effect of a damper on the 
wheel as follows. The angular velocity vector for the wheel, in a reference frame rotating with the 
wheel, has the components oj cos (X - o)t, co sin (X - a)t, and 12 + a. Comparing these with corre- 
sponding terms in the total angular velocity vector used in the derivation above (co cosXt, co sinXt, 12) 
shows that the wheel damper equations can be produced by changing the variables X and £2 such that 
X -*■ X - a and 12 -> £2 + a. Thus, 


and 


p'l + ~P\ + a,2 )'°i = (X _ S2 - 2 u)XqCo cos(X - a)t + ZqCO 2 

T' s = ot'xqCo[2p' 1 cos (X - o)t - p' x (X + 12) sin (X - a)t ] . 


( 6 . 6 ) 


In Equations 6.6 and in what follows, primed quantities are used for parameters associated with 
wheel dampers. Except for the difference in notation (p'j, Xq, and z' Q become z\ a\ and /j), the first of 
Equations 6.6 compares favorably with Equation 3.5, derived earlier. These equations can also be ana- 
lyzed to obtain an average value of T' s in a straightforward manner. 

For this particular damper orientation, then, we can find the average spin torque produced by the 
steady-state damper response to excitation arising from the general precessional motion of the system 
with an arbitrarily large cone angle. We now turn to other possible orientations of the dampers. 

F. Special Case 1 1 : Damper Normal to Spin Axis and Radius Vector 

In this case, let = e 3 = 0 and e 2 = 1 . Equations 6.3 and 6.4 then reduce to 


and 


c r k ~\ x o c ° 

p,+ — p.+ co 2 (l - sin 2 Xf) - 12 2 p , = - (12 + X)z A co sin X/ — sin2Xt 

1 m 1 m J 1 u 2 


(6.7) 


T s = m 


P i*o w 


pW 


- XqPj - 2p 1 p 1 12 + PjZqIX + 12)co cos Xt + — cos 2Xt + — - — sin 2 \t 


For a damper on the wheel, these equations become 


P'l + ^ Pi + ~ sin2 + CT ) 2 } P'i 

= - (12 + X)ZqCo sin (X - o)t - — — — sin 2(X - a)t 


*o “ 2 


and 


( 6 . 8 ) 


T' s = m' 


- XqPj - 2p' 1 p3(12 + ct) + P^z^X + 12)co cos(X - a)t 


p[x' 0 co 2 p' 2 co 2 

+ cos 2(X - a)t + — - — sin 2(X - a)t 
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G. Special Case III: Damper Radial 

In the case where e l = 1 and e 2 = e 3 = 0, the damper and torque equations become 

+ " cj2 ^ ” cos2 ^ - Pi = ~z 0 (X + £2)co cos Xt + x 0 (o> 2 sin 2 Xt + £2 2 ) 

2 


and 


T s = m 


P tO>‘ 


- 2p l (p l + x 0 )£2 - PjZgfX + £2)cu sinXt — ( 2x Q + p x ) sin 2Xt 


and 


For a damper on the wheel, these equations become 

p\ + ~'P'i + {^' -cj2[1 - cos 2 (X - CT)r] - (^2 + cr) 2 j p'j 

= “ZgCX + cos(X - a)t + Xg[co 2 sin 2 (X - a)t + (S2 + a) 2 ] 

T' s = m' - 2p' 1 (p' 1 + Xq)(£2 + a) - p\z' Q (X + £2)co sin (X - a)t 


(6.9) 


( 6 . 10 ) 


p' to 2 

— (2xg + p[) sin2(X- a)t . 

— 

The damper equations obtained for Special Cases II and III are clearly more complicated in gen- 
eral than those associated with Special Case I. They involve periodic coefficients and are similar in 
form to Mathieu equations with sinusoidal forcing functions of two frequencies. For small angles (co 2 
small), however, they reduce to tractable constant-coefficient equations, and they will be analyzed on 
that basis. 

The average spin-axis torque arising from the steady-state motion of the various dampers consid- 
ered above will now be determined, with the assumption of small-angle motion for the system. 

H. Small-Angle Analysis of the Special Cases 
Spacecraft Damper of Case I 

For Case I, the appropriate spacecraft damper equation of motion is 

c k 

p, + — p. + — p = (X - fi)x n co cos Xt , 

1 m 1 m 1 0 


with a steady-state solution 


p l = B' sin Xt + C' cos Xt . 
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The spin torque is given by 


i 

i' 


T s = mx Q co[2p 1 cos Xt - p t (X + £2) sinXt] , 
so that the average torque will be 


Since 


we have 


T s = mx 0 oj[B'X - VzB\\ + £2)] 
mx n co 

= —f~ 


0 

(x - n)x 0 

m 

k \ 2 

CO \ z 

m 


--x 2 




m 

m 



X- 

--X 2 



m 

m 



Vzcx qC 0 2 X(X - £2) 2 

(-^- x 2 ) 2 + ( x -^) 2 

\m / \ m/ 


Wheel Damper of Case I 

For Case I, the appropriate wheel damper equation of motion is 

Q t 

p\ + — p', + — , p, = (X - 12 - 2o)x n co cos(X - o)t , 

1 m 1 m 1 u 

with a steady-state solution 

p' x - B' sin (X - a)t + C' cos (X - a)t . 

The spin torque is given by 

T' s = m'xQCo[2pj cos (X - a)t - p'j(X + £2) sin (X - a)t ] , 


( 6 . 11 ) 
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so that the average torque will be 


T' s = m'x' 0 u[B'(\ -a)- ViB\X + £2)] 


Since 


we have 


m'x' n oj 

= — j-B\\-2o-Sl) . 


B' = 


0 


-Ck-o)—, 

m 


(X - £2 - 2a)x n oj ^-(X-a) 2 ! 
u m 


o) 2 -(X-a)-, 

m m 


(X - a) - 7 ~ (X - a) 2 

m m 


1 /2C , Xq 2 co 2 (X - a)(2a - X + £2) 2 


-(X-o) 2 

2 

+ 

(X - a) — 

m 


m 


Spacecraft Damper of Case II 

For Case II, the appropriate spacecraft damper equation of motion is 


P i 


+ — p, + (—- £2 2 ^ p, = -(£2 + X)z n co sin Xt , 
m 1 \m / 1 u 


with a steady-state solution 


( 6 . 12 ) 


Pj = B' sin Xt + C' cos Xt . 


The spin torque is given by 

T s = m[-x Q p 1 - 2p 1 p 1 £2 + p x z Q (X + £2)oo cosXt] , 
so that the average torque will be 


mz Q coC'(X + £2) 
2 
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Since 


C' = 


--ft 2 -X 2 -(fi + X)z n co 
m u 


— X 
m 


— - J2 2 - X 2 --X 


m 


m 


k 


we have 


T S = 


— x — - n 2 - x 2 

m 777 


1 /2Cz 2 co 2 X(X + n) 2 
\ra / \m / 


Wheel Damper of Case II 

For Case II, the appropriate wheel damper equation of motion is 


+ rn'P ' 1 


+ 


~ + o ) 2 

m 


p' x = - (£1 + X)ZqCO sin (X - a)t , 


with a steady-state solution 

/?j = B' sin (X - a)t + C' cos (X - a)t . 

The spin torque is given by 

T' s = m'[-x'Qp'i - 2p' 1 p' 1 (S2 + a) + p\z’ Q (\ + £2)co cos(X - a)t ] 
so that the average torque will be 

r s = V 2 m'z' 0 a>(\ + Q.)C' . 


Since 


C' = 


- (£1 + a) 2 - (X - a) 2 -(fi + X)z'co 
m u 


(X - o)— t 
m 


— r -(£l + o) 2 -(\-o ) 2 
m 


(X - a) , 
m 




A:' 


wz 


- - (£2 + a) 2 - (X - a) 2 


( 6 . 13 ) 
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we have 


1 /2C , z' ( ?U) 2 (\ - a)(£2 + X) 2 


n=r 


- (fl + a) 2 - (X - a) 2 

2 

+ 

(x-a)^7 

77?' 


m 


2 ‘ 


Spacecraft Damper of Case III 

For Case III, the appropriate spacecraft damper equation of motion is 

p, + — p, + (— ~ fi 2 )p, = ~zJX + £2)co cosXt +x n SX 2 
1 m 1 ' m / 1 u 

with a steady-state solution 

Pj = ,4' + 5' sin Xt + C cos X/ . 

The spin torque is given by 

T s = m[- 2p l (p l +x 0 )n - PjZ 0 (X + n)w sin Xt] , 
so that the average torque will be 


r c = - '/imz ft co(X + f2)S' . 


Since 


5' = 


X 

m 


~ z n (X + SX)co F2 2 - X 2 

0 77? 


^-fi 2 -X 2 X 

m m 


— x -~n 2 -x 2 

m m 


we have 




1 /2CZ 2 co 2 X(X + f2) 2 


(*_ n 2 -X 2 ) 2 + (-x ) 2 

\ 7?7 / \77? / 


Wheel Damper of Case III 

For Case III, the appropriate wheel damper equation of motion is 
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(6.14) 


(6.15) 



p[ + ~>p'\ + + o0 2 Jp'i =-z’qO. + ^)<^ cos (X - ff)f + Xq(S2 + a) 2 , 

with a steady-state solution 

p'.y = A' + B' sin (X - o)t + C cos (X - a)t . 

The spin torque is given by 

T' s = m'[~ 2 p[(p[ + Xq)(£2 + o) - p\z' 0 (\ + £2)co sin (X - cr)f ] , 
so that the average torque will be 

T' s = -V2m’z' 0 u(\ + Sl)B’ . 


Since 


S' = 


0 


~(X-a) 

m 


k' 


-z'(X + ft)cu ^-(ft + a^-CX-ay 
u m 


k' 


~ (fl + a) 2 - (X - a) 2 
m 


(X-a) 

m 


(X-a) 

m 


k! 


— : - (£2 + a) 2 - (X - a)" 
m 


we have 


1 /2m'z^w 2 (X-c;)(X + n) 2 


V = 
1 s 


7 ~ (£2 + a) 2 - (X- a) 2 


^7 (X-a) 
m 


2 * 


(6.16) 


I. General Observations 

Equations 6. 1 1 through 6. 1 6 are similar in many ways. Of particular interest is the preponder- 
ance of positive terms. Although the various equations differ in detail, it can easily be seen that the 
sign of X controls the sign of the average torque for each of the spacecraft dampers, and the sign of 
X-a determines the sign of the average torque for each of the wheel dampers. These quantities repre- 
sent the rotation rate of the transverse angular velocity vector in a reference frame fixed In the respec- 
tive body. In general, then, damping in a particular body increases the spacecraft spin rate if the trans- 
verse angular velocity vector moves counterclockwise in that body. For lack of a better word, this 
type of motion will be called a positive nutation rate. 

Earlier in this chapter, we showed that spin rate increases imply cone angle decreases. Thus, 
damping in a particular body decreases the cone angle if the nutation rate in that body is positive. For 
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a single-spin spacecraft, the nutation rate X, the spin rate £2, and the moments of inertia C and A are 
related by X = (C/A - 1)£2. The spin rate and the cone angle are related by Ci 1 = H cos 0, where H is 
the total angular momentum, which we take as a positive quantity that is a constant Of the motion for 
a system free of external torques. For 0 near 0°, stability would involve a decreasing cone angle, and 
for motion near 180°, an increasing 0 would indicate stability. In either case, stability can be seen to 
require C>A, the usual moment of inertia constraint for spinning bodies. 

Nutation (as defined above) is a measurable quantity in a spinning body, but it is not an obvious 
characteristic of torque-free motion to an outside observer. The spin rate is observable, of course, and 
if there is any cone angle present, another quantity, the precession rate, can be seen. This is the rate at 
which the nominal spin axis actually moves around the inertially fixed angular momentum vector, and 
it is one of the usual “Euler angle” rates in textbook descriptions of rigid body motion (\p 1 in Figure 
6.2). If the system angular momentum is expressed in terms of the Euler angle rate vectors, it can be 
easily shown that \p l = HI A, if 0 vanishes, as it does for the torque-free motion of a symmetric body. 

The relationship i p 1 = HI A holds for dual-spin systems as well, since the transverse angular mo- 
mentum is still Aoj. For the dual-spin system, we have the relationships 

CS2 + J' 3 o = H cos 0 , 


z 



Figure 6.2— Orientation of a rigid body by a 3-1-3 transformation of angles 

, d, and ip. 



and 




* l = 


// 


which can be combined to yield 


X = \p l cosd- SI . 

For the case of small cone angles, X ^ - SI. Spacecraft damping is stabilizing if X > 0, or in 
this case, if \f/ 1 > £ 2 ; wheel damping is stabilizing if X - a > 0, which in this case requires \p l > SI + o. 
Thus, we deduce that damping on a particular body is stabilizing if that body precesses faster than it 
spins. 

J. Energy Dissipation 

The damper equations just considered were all of the 

Pi + ~^l + K iPi =a + $i 
with a steady-state solution of the form 

Pj = A' + B' sin ( 7 1 + ipj) + C' cos ( 7 1 + < p^) . 

The energy dissipation rate for the damper is 

P = cp 2 , 


general form 
sin^r + v^) , 


and the average rate is given by 


P= l / 2 cy 2 (B' 2 + C' 2 ) . 


It can be easily shown that 


Thus 
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Term-by-term inspection of the average spin-axis torque expressions, Equations 6. 1 1 through 
6.16, shows that in every case they are of the form 



If a dual-spin system containing any arbitrary combination of the above types of dampers is con- 
sidered, all of the spacecraft dampers will be driven at the frequency X, and the wheel dampers at the 
frequency X- a. The average torques in the system will combine linearly, and we can say in general 
that 


p p 

ls XX -o’ 


(6.17) 


where P s and P w represent the rates of energy dissipation in the spacecraft and the wheel, respectively, 
and are inherently positive quantities. 

Since an increasing spin rate implies a decreasing cone angle, stability for an arbitrary initial state 
requires 

T s cos0q>O, (6.18) 

where 0 Q is either 0° or 180°. 

Equations 6.17 and 6.18 are of a form reminiscent of a result presented by Likins (Reference 13) 
after an “energy sink” analysis of the attitude stability of a similar dual-spin system. In fact, with 
allowances for differences in notation, the conclusions are identical. 

Working with linearized equations, Likins ignored equations “involving the changes in rotation 
rates of the two bodies about their common axis . . . under the assumption of small motor and bearing 
friction torques” and assumed that any damping present would decrease the kinetic energy of the sys- 
tem, i.e., T < 0. Then, recognizing the conservation of angular momentum in the system, he manipu- 
lated appropriate equations to deduce that the decreasing energy implies a stability criterion equivalent 
to that deduced here. 

Now, with the assumption of a constant wheel rate, it can be shown that in the case of prime 
interest in dual-spin stabilization (i.e., the stability of pure rotations about the axis of minimum mo- 
ment of inertia), stable equilibrium involves a maximum energy state. That is, stable behavior implies 
an increasing kinetic energy for the system. This can be seen by combination of the equations 

T = YiAco 2 + 'MC-J’^n 2 + VzJ'^n + a) 2 , 


and 


Aco = H sin 9 , 
CYl +/3U -H cos 9 
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to obtain 


T = ^AC Sin2 ® + constant * 


(6.19) 


Equation 6.19 shows that, for the dual-spin system, 0°, 90°, and 180° will always represent ex- 
trema in kinetic energy. For “favorable” systems (i.e., C>A) pure spins about the nominal spin axis 
will be minimum-energy states. In this, the dual-spin system is similar to a single-spin body, but sta- 
ui^ equilibrium might involve any cone angle (0° < < 180°). Thus, in terms of kinetic energy, max- 

imum, minimum, or even intermediate energy states can be stable. 

In light of the above considerations, the previously mentioned agreement with Likins’ result is 
somewhat surprising. In fact, one might reasonably reject a priori “energy sink” analyses in the area 
of dual-spin stabilization. 

However, it has been shown that our stability criterion takes on a simple form when expressed in 
terms of the amount of energy dissipated on each of the bodies, although we have not said anything 
about the overall system kinetic energy. Our result, as given by Equation 6.17, can be interpreted in 
terms of a dimensionless stability diagram [in the (o/£l 0 )q -plane representation used earlier]. 

Again, for definiteness of cone angle only, it is assumed that £2 Q > 0. The initial-state plane (Fig- 
ure 6.3) then becomes divided into various regions characterized by the signs of the quantities cos0 Q , 


Q 



59 



X, and X - a. Quantities pertinent to the subdivision of the initial-state plane shown are repeated here 
for reference. 


Q 


C 
A ’ 



X = (q - l)fi 0 + q'o , 

and 

H 

q£l 0 + <?'a = -cos0 o . 

With the previously mentioned assumption that > 0, these regions represent the combination 
of variables shown in Table 6.1 . 

Regions 3 and 4 are seen to be always stable and always unstable, respectively. In the other re- 
gions, however, stability or instability depends on the relative dissipation in the two bodies, and the 
requirement for stability becomes 

P s X 

(1) — > r for Region 1. 

P W °~ X 

p s X 

(2) — < for Region 2. 

P W °~ X 

p s X 

(3) — > — — r for Region 3. 

P pj / Q “ 

Thus, the quantity k x = X/(a - X) represents either an upper or lower limit for the ratio P s /P w for 
stability. Lines of constant k l can be determined algebraically and simply as follows: 


Table 6.1— Characteristics of initial state regions. 






Effect of Damping 

Region 

Cone Angle 

Sgn (X) 

Sgn (X - a) 



Spacecraft 

Wheel 

i 

180° 

— 

+ 

Stabilizing 

Destabilizing 

2 

0° 

— 

+ 

Destabilizing 

Stabilizing 

3 

0° 

+ 

+ 

Stabilizing 

Stabilizing 

4 

0° 

— 

— 

Destabilizing 

Destabilizing 

5 

0° 

+ 

— 

Stabilizing 

Destabilizing 
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but 


Thus 


or 


X , ^ . *i 

= k, => X = , . , a 


ff-X 1 


1 +k. 


X = (q-l)n 0 + q'a 




1 + ^1^0 


<7 = 1 + 



These lines all pass through the point (0, 1) and are obviously straight lines with a slope of 

VO +VV- 

In the limiting cases of = 0 and k l = «>, these lines are q + q’(o/£l Q ) = 1 and q + q'(o/£l Q ) = 

1 + o/£l 0 . The line will be horizontal when k 1 = q'/(l - q'). Figure 6.4 shows this overall result. 

The above result suggests that the transformation of the stability diagram, discussed earlier in 
terms of the parameters 5 and 5' for the specific system analyzed, can be understood quite simply in 
terms of the quantity k 1 , which measures the relative energy dissipation in the spacecraft and the 
wheel. Figure 6.5 shows how the transformation progresses. 
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Figure 6.5— Transformation of stability diagram in terms of k^. 
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CHAPTER 7 


SUMMARY AND CONCLUDING REMARKS 

The research reported here concerns the attitude stability of a dual-spin spacecraft whose behav- 
ior is described by a seventh-order, nonlinear, nonautonomous set of ordinary differential equations, 
inconveniently coupled even in the highest derivatives. 

In the recent literature, this and similar systems are analyzed by linearization of the equations of 
motion about a reference motion that is a pure rotation about the nominal spin axis, or equivalently, 
linearization in six of the seven dependent variables. One of the resulting equations integrates trivially, 
reducing the order to six. The remaining equations have periodic coefficients and are studied by nu- 
merical Floquet analyses of very special cases. 

With damping on one body only, however, the linearized equations have constant coefficients, 
and the order of the system is reduced to four. More general results are then obtainable by a Routhian 
analysis of the characteristic equation for the system. Five Routhian inequalities necessary and suffi- 
cient for stability are produced by this procedure. Two are trivially satisfied for practical systems, one 
is a somewhat general criterion not related to the damper mass sizes, and the other two impose damper 
size limitations. 

The initial motivation for this work was a curiosity about the nature of instabilities in such a 
dual-spin system. Unstable single-spin spacecraft are subject to the well-known flat-spin or tumble in- 
stability, but nothing could be found in the literature concerning analogous behavior for dual-spin 
systems. 

In the approach to this problem, gross linearization of the equations of motion was rejected 
a priori , since “large angle” motion of the system was anticipated where instabilities occur. Physical 
insight eventually led to the perturbation technique applied to this problem. 

Through consideration of the nonlinear equations of motion with the damper masses vanishingly 
small, the “rigid body” behavior of the system was revealed, and the unperturbed motion of the sys- 
tem was found to be similar to the general precession of the torque-free, symmetric, rigid body in 
classical mechanics. Then, with the retention of all terms up to the first order in the damper masses, 
equations were produced that represented the first-order behavior of the dampers, and these were 
solvable in closed form. 

Since an inertially fixed angular momentum vector is necessarily a constant of the motion for the 
system, in spite of the complicated nature of the equations, the spin component of this momentum 
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and the equations of motion were used to deduce an exact differential equation controlling the behav- 
ior of the system cone angle. This equation is of a form amenable to analysis by an averaging tech- 
nique, which produced an equation of the form 6 =f(6) governing secular changes in the. cone angle. 

Since in spin-stabilized systems the behavior of the cone angle alone often represents all one 
wishes to know about stability, the perturbation analysis had, in effect, determined the essential be- 
havior of a system characterized by a seventh-order, nonlinear, nonautonomous set of differential 
equations with arbitrary initial conditions in terms of a first-order equation. Plotting the function f(d) 
for a particular configuration immediately reveals the stability characteristics of the system, indicates 
“zones of attraction” for particular stable equilibrium points, and shows the “time constant” associ- 
ated with cone-angle buildup or decay near equilibrium states. 

In appropriately reduced form, the general stability criterion developed was compared to perti- 
nent results in the recent literature, all of which treat “small-angle” stability as mentioned before, and 
generally good agreement was found. In addition, however, analytical stability boundaries were estab- 
lished where only numerical Floquet results were found in the literature, and the nature of possible 
instabilities was revealed. 

Two general types of instability were found-complete inversions of the system and an eventual 
steady-state cone angle. In the former case, the spacecraft literally turns upside down in reaching a 
stable equilibrium state, and in the latter, calculable stable cone angles ranging from 0° to 180° are 
found. The inversion phenomenen is unreported in the literature, if known at all, and the flat-spin 
instability of single-spin spacecraft is a special case of the other type of instability. 

Fairly general stability diagrams appear in the literature for the special cases of damping on only 
one of the bodies. These diagrams are shown to be consistent with the general first-order stability 
criteria developed here, and, in addition, unstable regions are divided according to the nature of the 
instability. Also, these diagrams are generalized to cases with damping on both bodies, and the trans- 
formation of the stability diagram as damping is gradually shifted from one of the bodies to the other 
is revealed. 

Certain possible instabilities revealed by Routhian and Floquet analyses in the literature were 
undetected in the perturbation analysis reported here. These are characterized by inequalities which 
limit the sizes of the damper masses. The significance of these limitations is greatly enhanced in cases 
where one attempts to “tune” the dampers to very low excitation frequencies. In practice, these in- 
stabilities are unlikely, due to the unrealistically low spring constants required for tuning. Since the 
pertinent inequalities are related to the sizes of the damper masses, it is felt that they represent second- 
order effects relative to the perturbation analysis reported here. They arise because of the coupling 
between equations and were undetected here because the equations uncoupled in this firshorder 
perturbation analysis. 

Possible generalizations of the results obtained are sought by a determination of the effect of an 
arbitrarily positioned and arbitrarily oriented spring-mass-dashpot damper in a dual-spin system. In 
general, the damper equations derived are of the Mathieu type, so that instabilities within the dampers 
themselves can occur. For small-angle motion, however, all of the equations reduce to constant- 
coefficient linear equations solvable in closed form, and the spin-axis torques generated by their 
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motions can be found. Although quantitative differences exist for various damper locations and orien- 
tations, all have in common the characteristic that damping on a particular body is stabilizing if that 
body precesses faster than it spins. For a single-spin body, this statement is equivalent to the maximum- 
inertia axis rule for stability of rotations in the presence of energy dissipation. 

The average spin-axis torque generated by a particular damper is found to be simply related to 
the average rate of energy dissipation in that damper and the excitation rate at which it is driven. This 
fact permits a generalization to a simple stability criterion applicable to a system with any number of 
arbitrarily positioned and arbitrarily oriented spring mass-dashpot dampers. This criterion is consist- 
ent with a result obtained in the literature by an “energy sink” approach to dual-spin stabilization, but 
the agreement is, in principle, puzzling. The energy sink approach assumes a monotonically decreasing 
kinetic energy, but, in the case of our system, stable behavior often involves a move to a maximum 
energy state. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland, May 29, 1970 
877-1 1-75-01-5 1 
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Appendix I 


Derivation of the Equations of Motion for the System 


First, the Lagrangian will be determined for the system. It will then be substituted into Lagrange’s 
equations of motion with the angular velocities in the spacecraft-fixed coordinate system as quasi- 
generalized velocities and the displacements of the dampers’ active masses as generalized coordinates: 


and 


d ( d£\ d£ d£ 

+ CO. co. 


dt\duJ ' 2 9co 3 3 9co 2 


d ( b£\ _d£_ _ d£ 

dt \9co 2 ) °° 3 9cOj C ° 1 9co 3 


d_ I b£ \ 9£ _ d£ 

dt \9coJ Wl 9co 2 L ° 2 9co 3 


dt 9£\_9£ 
dt \9z / dz z ’ 

d_ /9£\ _ d£ = 
dt \9z 7 9z' Qz ' - 


= 0 , 
= 0 , 
= 0 , 


The kinetic energy T will be developed in three steps: 

T 0 = kinetic energy of the system with z = i = z' = z' = 0, rotating about and translating with 
respect to the true center of mass. 

T 1 = additional kinetic energy of the spacecraft damper active mass due to z and z ¥= 0. 

T 2 = additional kinetic energy of the wheel damper active mass due to z' and z' =£ 0. 

The total kinetic energy will then be T = T Q + + T v 

The center of mass of the system will lie along the nominal spin axis at a distance z above the 
spacecraft center of mass, as depicted in Figure 1. 1 . 

Since the total mass of the wheel is M' + 4 m' b and the total mass of the spacecraft is M + 4 m b ,z 
becomes 
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Figure 1.1— System center of mass. 


(M' + 4 m* b - m f )l + m'(7 + z') + mz 


Z " A/ r 

where 

A/^ = A/ + 4w 6 + M' + 4m’ b 

is the total mass of the system. 


If 

M + 4m, 

L — ~ ~ / 

1 Af r 

and 

AT + 4m' 
'3= *'■ 

then 

mz + m'z* 

z = /, + 

2 m t 


= / 2 + ^ 

and 

1 

N 1 
II 

— i 

1 




where 


m 



and 

? = pz + p'z' . 

Translation of the system relative to the center of mass is along the nominal spin axis at a rate f . 
Thus, the translational kinetic energy is simply Vi M T $ 2 . 

The transverse moment of inertia of the nominal system about the true center of mass is 
/ x + 2 m b a 2 + (M + 4 m h )(l 2 + f) 2 + /' + 2 m' b a t2 + (M* + 4 m' b )(l x -$) 2 =A + M T $ 2 . 

In the above equation, 

A =/ 1 +/j + 2m b a 2 + 2 m' b a t2 + (M' + 4m^)(l - r^)/ 2 , 

where 


M' + 4m; 



The transverse kinetic energy of rotation of the system is thus Vi(A + M T $ 2 )(co 2 + co^)- 
The kinetic energy of rotation about the nominal spin axis is 

Vi (/ 3 + 4m b a 2 )Gj 2 + l A(I\ + 4m;a ,2 )(cc 3 + a) 2 = ViCco 2 + + a) , 

where 

C = / 3 + 4m b a 2 + / 3 + 4m;<s' 2 

and 

= /; + 4m;a' 2 . 

Combining the above quantities yields 

r 0 = V 2 MJ 2 + y 2 Cco 2 + 1 / 2 /;a(2co 3 + a) + %(i4 +M r f 2 )(co 2 + cu 2 ) . 

Determination of T x 

Now, Tj = Vim(v 2 - Vq), where v is the actual velocity of the active mass of the spacecraft damper 
and v 0 is the velocity of the same mass if z = z = 0. Let v m be the position vector for the mass m. Ex- 
pressed as a column matrix, 
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Figure 1.2 shows the positions of the mass m, the system center of mass, and the total angular 
velocity vector oj t in a spacecraft-fixed coordinate frame. 

Now, 


Thus, 


and 


Then, 


where 


03 T* r m = 


W 1 


CO, 


a 0 z - L - $ 



v = v n + 5v , 



5v = I -zco 


and 


Hence, 

v 2 - Vq = 2v 0 ■ 5v + (5v) 2 

= 2{[co 2 (-/ 2 -f)zco 2 ] + [a to 2 -tOj(-/ 2 - flK-ztOj) + (-f)z } + z 2 co 2 + z 2 co 2 + z 2 
T 1 = i /2ot|2co 2 z[-co 2 (/ 2 + f)] + co^z 2 - 2co 1 z[flto 3 + + f)] + co 2 z 2 - 2z(f +aco 2 ) + z 2 | . 
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Determination of T 2 

Now, T 2 = x /i m'(v' 2 - v' Q 2 ), where v' is the actual velocity of the active mass of the wheel damper 
and Vq is the velocity of the same mass if z' = z' = 0. Let r' be the position vector for the mass m' 
(Figure 1.3). Then, 



and 


o> r x r 


1 

co, 


co- 


co. 


a' cos t jj a' sin \p z' + 1 j - f 


Thus, 


v = 


vk - 


/-a' a sin <// + co 2 (z' + /j - f) - c o 3 a' sin t// N 
a' a cos \p + co 3 a' cos \p - co 3 (z' + l x - f) J , 
z' - £ + c o^a' sin - c o 2 a’ cos \p 

■a'a sin i p + co 2 (/ 1 - f) - c o 3 a' sin \p^ 
a' a cos i p + co 3 a' cos \jj - co t (/j - f) J , 

-f + co 3 a' sin \p - co 2 a' cos \p 
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and 


v' = Vq + Sv' , 


where 

Thus, 

v' 2 - v' Q 2 = 2vq • 5v' + (5v') 2 

= 2|[-a'a sin 1// + w 2 ^i ” £) - 03 $ a ' s ^ n ^](z'co 2 ) 

+ [a'o cos \p + co 3 a' cos \ p - coj (/j - JOK-z'tOj) 

+ (- j + cojfl' sin \l> - c o 2 a' cos i//)z'j + z' 2 co 2 + z' 2 co 2 + z' 2 . 

Thus, T 2 becomes 

T 2 = }4ra' |2z'co 2 [-a'a sin \p + co^/j " f) “ sin + z' 2 co 2 
- 2z'co 1 [a'a cos \p + co 3 a' cos \p - cojC/j - f)] + z' 2 co 2 
+ 2z'[~i + cojfl' sin i// - to 2 a' cos i//] + z' 2 | , 

and the total kinetic energy of the system is 

T = ViM T t 2 + V 2 C 0 J 2 + y 2 J' 3 o(2u 3 + o) + V 2 (A -M T t 2 )(oj 2 + co 2 ) 

+ 1 /277j|2C0 2 z[-C0 2 (/ 2 + f)] + <^ 2 Z 2 - 2c0jZ[flC0 3 + COj(/ 2 + $■)] + C0 2 Z 
- 2z(t + aco 2 ) + z 2 j + 1 /2m'|2z'co 2 [-fl'a sin 1 p + co 2 (/j -f) - co 3 a' sin \ p] 

+ z' 2 co 2 - 2z'co 1 [a'a cos \p + c o 3 a' cos 1// - C0j(/j - £)] + z' 2 co 2 
+ 2z'( - f + cOjfl' sin 1/7 _ co 2 a' cos + z' 2 | . 

The potential energy of the system is that energy stored in the damper springs, i.e., V = Vi^z 2 + 
y 2 k'z' 2 ; and the generalized forces associated with the coordinates z and z' are the viscous drag forces 
Q z = -cz and Q z , = -c'z'. The Lagrangian is, of course, £ = T + V, and the partial derivatives (and their 
time derivatives) necessary for construction of the equations of motion for the system are as follows: 
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d£ _ dT 
3 gj 1 doo l 

= (A + M t £ 2 )gj x + mjcjjlz 2 - 2 z(l 2 + f)] -zaco 3 j - m'[z f a ' cos \jj (co 3 + a) - z'a! sin \ p] 

+ m'o l [z ,2 + 2z'(l l -?)] , 

7 pv n 

JidZT = \ A + M ^ 2 + m[z2 ~ 27(/ 2 + m + m ' [z ' 2 + 2z ' (/ i " £)]}*! + coj {2 m t k 

+ m[2zz - 2 zf - 2 z(/ 2 + f)] + /n'[2z'z' - 2z'f + 2i'(/j - f)]} + m'a' sin \}/ z' - mazai^ 

- mazoj^ - m'lz'a' cos to 3 - z'a'o cos i// + (co 3 H* a)(z V cos \p - z'a’o sin i//)] , 

ax _ ar 

3co 2 3o> 2 

= [z4 + M r f 2 + mz 2 - 2 mz(/ 2 + f) + m'z' 2 + 2m'z'(/j - f)]a; 2 ” “ tn'a'z' sin \p (co 3 + a) 

-m f a f cos \p z' , 


d 3jC 
dt 3co 2 


= [y4 + A/ r f 2 + mz 2 - 2 mz(J 2 + f ) + m'z' 2 + 2 m'z'(/ 1 - f)]co 2 + co 2 [2Af r ff + 2mzz 

- 2mzf - 2 mz(/ 2 + f) + 2m'z'z' + 2m'i'(/ 1 - f) - 2 m'z'f] - m<zz - mVz' sin t// co 3 

- mV(co 3 + a)(z' sin i// + z'a cos i//) - m'a' cos \p z f + m'a'z'a sin i// , 


3jC = 3T 
3cc 3 3co 3 

= Cco 3 + / 3 a - mazco l - m'a'z'( co 2 sin i// + co t cos \I/) , 
d 3jC 

— = Cco 3 - ma(zd) l + zco x ) - m Vz'( co 2 sin ^ + co t cos i//) - m Vz'(sin \p 6 j 2 - cos i// aoo 2 

+ cos \p cb 1 - sin i// acoj) , 


d£ _ 3T 
3z 3i 

= -m'z'p + mi - mpz - maco 2 , 


d_ d£ 
dt 3 z 


= “m'z'p + m(l - p)z - maco 2 
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d£ _ dT_ _ dV 
d z dz d z 

= mz(co 2 + to 2 ) - m(l 2 + p'z'Xce^ + co|) - mpz((x>\ + oo 2 ) - maco { co 3 - kz , 

9 £ _ 97 ; 
dz' dz' 

- m'z' + m'(- pz - p'z' + sin i// - co 2 a' cos i//) , 
d 

— — = m'[(l -p')z' - pz + cc.a' sin < // - co.a' cos \p a - co 9 a' cos \p + cocci' sin i p a] , 
at oz 11 z 

d£^dT_dV 
dz' dz' dz' 

= m'tc 2 [-a'a sin i// + co 2 (/j - pz) - c o 3 a' sin \jj] - m’p’z'oa\ + m'z'(co 2 + to 2 ) - m'p'z' oj 2 
- m'ccja'a cos \p + ce 3 a' cos \p - co^/j - pz)] . 

The resulting equations of motion are 

[. A +M t ^ 2 + mz 2 - 2m(? + l 2 )z + m'z' 2 - 2m'(? - / J )z']co 1 - (maz + m'a'z' cos \p) co 3 + m'a' 
= 04 -C)cu 2 co 3 -J' 3 au> 2 -2M t ^uj 1 +A/ r f 2 co 2 co 3 + 2m(? + / 2 XZC 0 J -zco 2 cc 3 ) 

+ mz(2£aJ 1 - 2 ic 0 j + zco 2 cj 3 + atOjte 2 ) + 2m'(? “ / 1 )(z'co 1 -z'co 2 <o 3 ) 

+ m'z'(2£co x - 2z'co 1 + z'cj 2 oj 3 + a' cos \p cOjCJj) - m'a' sin i// z'[(to 3 + ct) 2 - co 2 ] , 

[^4 + Myf 2 + mz 2 - 2m(f + l 2 )z + m'z' 2 - 2m' - / 1 )z']co 2 - m'a' sin <// z'co 3 - muz - m'a' cc 
= (C -/liojjOjg +J' 3 au> x - 2M T ^u> 2 -M T $ 2 co 1 co 3 + 2m(? + / 2 )(zct> 2 + zt 0 jCJ 3 ) 

+ mz(2fa> 2 - 2 zco 2 - ztOjO^) + muz(co 2 - to 2 ) + 2m'(J - / 1 )(z'cj 1 cj 3 + z'co 2 ) 

+ m'z'( 2fcJ 2 - 2z'co 2 -z'cOjCo 3 ) + m'z'[(co 3 + a) 2 - to 2 ]a' cos i// - m'a' sin \]j z'co 1 co 2 , 

-( maz + m'a' z' cos i//)co, - m'a' sin t // z'to 2 + Cco 3 

= ma(2zco, - z oo 2 a> 3 ) + m'a' sin \p (2z'co 2 + z'cojtOj) + m'a' cos i // (2z'co 1 - z'co 2 co 3 ) , 

-mate 2 +m(l -p)z -m'pz' = -matOjtOj + m(te 2 + to 2 )[z(l -p) ~/ 2 - p'z'] -cz -kz , 


- k'z' 

sin \p z' 

( 2 . 1 ) 
is i// z' 

( 2 . 2 ) 

( 2 . 3 ) 

( 2 . 4 ) 
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and 


m'a' sin \ jj cb 1 - m’a' cos \jj 6 j 2 - mp'z + m'( 1 - p')z' 


- -m'a'[sin i p co 2 (co 3 + 2a) + cos \ p gJj(co 3 + 2a)] 
+ m '{ coj + oo 2 )[z'(l ~p') + /j -pz] - c'z' - k'z' . 


( 2 . 5 ) 
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Appendix II 


Damper Size Limitations Imposed by Routhian 
Analyses When Dampers Are Tuned 


In Reference 1 1, Mingori explored a parameter space with 5 = ma 2 lA and 5' = m f a' 2 /A in the 
range from 0.001 to 0.100. Of additional interest, however, is the behavior of the system when one or 
the other of these parameters vanishes, i.e., on the 5 and 5' axes. The results of Mingori’s Routhian 
analysis are then applicable, and we deduce that three of the five inequalities produced by this analysis 
are not of particular concern to us here. Two of them are trivially satisfied for practical systems, and 
the other represents the stability criteria which have been shown to be in agreement with that deter- 
mined for these special cases in the research reported here. 

The following inequalities, however, must also be satisfied for stability, and they depend on the 
size of the damper masses. In Mingori’s notation, these are 

K+ 5 2 (A 2 -r£-Ar n )> 0 (II. 1) 

and 

A(AA - 5 2 4) > 0 , (II. 2) 

for damping on the spacecraft only, and 

K’ + 5' 2 [(A - rf - (A + r n )(r CT + r n )] > 0 (II.3) 

and 

(A-/- a )[(A-r a )r-5' 2 (r n + r a ) 3 ]>0, (II.4) 


for damping on the wheel only. 

Before attempting a Floquet analysis of the general case (i.e., damping on both bodies), Mingori 
made the assumptions that p = 5 and p' = 5' and tuned the dampers such that 

A - = 8(1 -5)A 2 (II. 5) 


and 


K’ = 5'(1 -8')(A-r o ) 2 . 


(II.6) 


Substituting Equations II.5 and II. 6 into Equations II. 1 through II.4 produces the modified 
Routhian inequalities 
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A 2 - 8r n (r n + A) > 0 , 

(II.7) 


A 4 - 5A(A 3 + 4) > 0 , 

(H.8) 


(A - rf - 5'(A + r n )(r a + r n )> 0 , 

(II.9) 

and 

(A - rf - 5'[(A - rf + (r n + rf ] > 0 , 

(II. 10) 

where 

(q - \ )£l 0 + q'o 

A = 

n 






v ’ 



>! 

II 



_ a 


and 

r °~ r? ’ 

V = l<7^ 0 + • 


In terms of more familiar quantities, we have 



X 3 - 5fi 0 (n 0 + X 0 ) > 0 , 

(11.11) 


A 4 -5A 0 (A 3 +ft 3 )>0, 

(11.12) 


(X 0 - of - 5'(A 0 + n Q )(n 0 + a) > 0 , 

(11.13) 

and 

(A 0 - a) 4 - 5'(A 0 - a)[(A 0 - of + (ty, + a) 3 ] > 0 . 

(11.14) 


Rearranging to see the corresponding damper size limitations, we have 
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Appendix III 


Derivation of Quadratic When F = 0 


Consider. 


Let 


Then 


X 0 (2o - X 0 + S2 0 ) = ±(3(o- X 0 )(X - S2 0 ) . 
x 0 = (q ~ 1 )S2 0 + q'o . 




1 ]K + *" 2 ) 


or 


[q'(2 - q') ± (3q'(q' - 1 )] (— ) + [(q - 1 )(2 - q') - q\q - 2) ± (3q'(q - 1 ) 


± 00 -q')(2-q)] T - + (q-\)(2-q)±m -q)(2~q) = 0. 
“o 


For q' = 0.2, tliis becomes 


(0.36 + 0. 1 6/3) + [ 1 .6? - 1 .4 ± 0( 1 .4 - 0.6 q)] + ( 1 ± (3)(q ~l)(2~q) = 0 
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Appendix IV 


Derivation of Equation 6.4 

F x p - r Q x (F - F q ) = mj- 2 p x (cj r x p) - X(r • k)p x cj + X(r • ai )p x k 

- (oJ r • r)p x co r + oj^p x r - r Q x p - 2r Q x (co r x p) 

- X(p • k)r Q x a) + X(p • a>)r 0 x k - (oi T ■ p)r Q x o? r + cj^.r 0 x pj 
= m|-2p 1 p 1 co r + 2 (p • CJ r )p - X(r • k)p x to + X(r • co)p x k 

- (co T ■ r)p xw r + w^p x r - r 0 x p - 2(r 0 • p)o> r + 2(r Q • co T )p 

- X(p • k)r 0 x cj + X(p • co)r Q x k - (cc T • p)r Q x oj t + oj^r Q x p| . 

k • [F xp + r Q x (F-F 0 )] = m|-2p 1 p 1 fi + 2 (p • co r )Pje 3 -PjX^ + p 1 e 3 )0 1 GJ 2 -e 2 Wj) 

- (co r • rXejCOj -e 2 oj 1 )p 1 -e 2 x Q p x -p 1 Xe 3 x 0 co 2 

— (ojy ■ p)Xq(x >2 ~ 2(Yq • p)£2 + 2(r 0 • c«Jj.)6 3 p 3 J . 
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